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Zero-field muon spin relaxation experiments have been carried out in the Pr(0si_2,Ru2:)4Sbi2 
and Pri_j;LajiOs4Sbi2 alloy systems to investigate broken time-reversal symmetry (TRS) in the 
superconducting state, signaled by the onset of a spontaneous static local magnetic field Bs . In both 
alloy series Bs initially decreases linearly with solute concentration. Ru doping is considerably more 
efficient than La doping, with a ~50% faster initial decrease. The data suggest that broken TRS is 
suppressed for Ru concentration x > 0.6, but persists for essentially all La concentrations. Our data 
support a crystal-field excitonic Cooper pairing mechanism for TRS-breaking superconductivity. 

PACS numbers: 71.27.+a 74.70. Tx 74.25.Nf 75.30.Mb 76.75. 
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In unconventional superconductors symmetries in ad- 
dition to gauge symmetry are broken in the supercon- 
ducting state, leading to novel properties and the possi- 
bility of more than one superconducting phase. -'^ Breaking 
of time-reversal symmetry (TRS) by a superconducting 
transition is a relatively rare example of such additional 
broken symmetry. Strong experimental evidence for bro- 
ken TRS comes from zero-field muon spin relaxation (ZF- 
/iSR) experiments that observe the onset of a sponta- 
neous local field Bs below the superconducting transition 
temperature Tc- Spontaneous fields have been observed 
by ZF-^SR in (U,Th)Bei3,2 UPts^ (although not with- 
out controversy"*'^), Sr2Ru04,^ the first Pr-based heavy- 
fermion superconductor PrOs4Sbi2)^ LaNiC2^ and, re- 
cently, PrPt4Gei2.- The ZF-^SR techniquCfiS in which 
spin-polarized muons are stopped in the sample and pre- 
cess in their local fields, is very sensitive to small static 
fields and thus is ideally suited for the study of broken 
TRS in superconductors. 

The isostructural fiUed-skutterudite compounds 
PrRu4Sbi2 and LaOs4Sbi2 are both conventional 
BCS-fike superconductors (T^ = 1.1 K and 0.74 K, 
respectively) f^ii^ Superconductivity is found for all 
values of Ru or La concentration in the alloy series 
Pr(Osi_a;Rua;)4Sbi2^ and Pri_yLayOs4Sbi2,^4 with 
relatively slow changes of Tc with composition. This 
is quite different from the behavior of the majority of 
heavy-fcrmion superconductors, where chemical substi- 
tution rapidly suppresses Tc- In Pr(Osi_a;Rua;)4Sbi2 Tc 
decreases smoothly from 1.85 K at a; = to a minimum 



of ^0.75 K at a; w 0.6, and then increases to 1.1 K 
at a; = 1 (Ref. [isl ). In Pri_j,Laj,Os4Sbi2 Tc decreases 
monotonically with y to 0.74 K at y = 1 (RcL 14). This 
behavior raises the question of how the TRS-breaking 
superconductivity of PrOs4Sbi2 evolves with Ru and La 
substitution. 

This Letter reports the results of ZF-/iSR exper- 
iments in Pr(Osi_a;Rua;)4Sbi2 and Pri_j,Laj,Os4Sbi2, 
which were undertaken to study the evolution of Bg with 
Ru and La doping. Preliminary results have been re- 
ported previously.*'^ An initial linear decrease of Bs with 
solute concentration is observed for both alloy series, but 
the data suggest very different effects of Ru and La: Bs 
is suppressed ~50% faster by Ru doping than La dop- 
ing and extrapolates to zero near the minimum in Tc{x) 
(Ref. d3i) , whereas for La doping broken TRS appears to 
be present for most if not all La concentrations. Our 
results support the theory of TRS-breaking supercon- 
ductivity from pairing via itinerant crystal-field cxcita- 
tionS )*^d'' and motivate further studies of these systems. 

The samples of PrOs4Sbi2, Pr(Osi_^Ruj;)4Sbi2, and 
Pri_yLaj,Os4Sbi2 used in this study consist of randomly- 
oriented small (^0.1 mm) crystallites prepared by the 
Sb-flux method. Strong de Haas-van Alphen signals ob- 
tained from similarly-prepared crystals** attest to their 
high quality. ZF-/iSR experiments were carried out at 
the Meson Science Laboratory, KEK, Tsukuba, Japan, 
and at the ISIS Neutron and Muon Facility, Rutherford 
Appleton Laboratory, Chilton, U.K. 

Figure[l]shows the time evolution of the decay positron 



2 



count rate asymmetry, proportional to the positive-muon 
(^+) spin polarization (Ref. [13), in PrOs4Sbi2 

and representative alloys at temperatures above and be- 
low Tc- A constant background signal originating from 



0.2 




Time (iis) 

FIG. 1. (Color online) Time evolution of muon decay positron 
asymmetry, proportional to the muon spin polarization P^{t), 
above and below the superconducting transition in PrOs4Sbi2 
and representative Ru- and La-based alloys. 



muons stopping in the sample holder has been subtracted 
from the data. As previously reported,^ in the end com- 
pound PrOs4Sbi2 the relaxation becomes faster in the 
superconducting state. Similar increases are observed in 
the alloys. 

The ZF-/iSR data are well described by the damped 
Gaussian Kubo-Toyabe (K-T) functioniii^ 

P^(i) = cxp(-At)GfT(A,<), (1) 

where 

^(A,i) = J + ^(1 - A2i')exp(-iA2i2) (2) 

is the K-T functional form expected from an isotropic 
Gaussian distribution of randomly-oriented static (or 
quasistatic) local fields at muon sites. The rms width 
of the static field distribution is where 7^ is the 

muon gyromagnetic ratio, and A is the exponential re- 
laxation rate associated with an additional contribution 
to the muon spin relaxation. Both A and A contribute 
to the increased low-temperature relaxation (Fig. [T]). 



Figure [2] shows the temperature dependence of A in 
Pr(Osi_:,Ru:r)4Sbi2, a; = 0, 0.1, 0.2, and 0.3, and 
Pri_j^Laj;Os4Sbi2, y = 0, 0.2, 0.4, 0.6. In the normal 
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FIG. 2. (Color online) Points: temperature depen- 
dence of the ZF Kubo-Toyabe static relaxation rate A in 
(a) Pr(0si_iRu2:)4Sbi2 and (b) Pri_j;LayOs4Sbi2. Curves: 
fits of Eq. ((3]) to the data using the temperature dependence 
of the BCS order parameter for Ae(T') (barely visible for 
x,y = 0). Arrows: Tc from bulk measurements. 



state A is due to dipolar fields from neighboring nuclear 
magnetic moments. An increase in A below T,. is ob- 
served in in both alloy series, indicating the onset of a 
spontaneous field in the superconducting state. The size 
of the increase becomes smaller with increasing solute 
concentration. To within errors no increase is observed in 
the end compounds PrRu4Sbi2^^'^° and LaOs4Sbi2'^'^ 
In superconducting PrOs4Sbi2i and Pro.8Lao.20s4Sbi2 
longitudinal applied fields > Hd (50-100 Oe) "decou- 
ple" the K-T relaxation, indicating that it is indeed qua- 
sistaticpiS 

Below Tc the nuclear dipolar and electronic contribu- 
tions to A are uncorrelated and add in quadraturei^ 

A{T) = [Al + Al{T)Y'\ (3) 

where A„ is the normal-state nuclear dipolar rate and 
Ae(T) is the additional relaxation rate due to the spon- 
taneous field from superconducting electrons. The K-T 
form assumes this field, like the nuclear dipolar field, is 
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randomly distributed, although as previously noted^ the 
data cannot discriminate between random and uniform 
spontaneous fields. A spontaneous internal field from 
broken TRS is expected only if (a) the superconductor is 
inhomogeneous and the field is nonuniform,^ or (b) the 
pairing is nonunitary and the probe spin is hyperfine- 
coupled to the (uniform) Cooper-pair spinj^ Thus it is 
not possible to decide between these alternatives from 
the ZF-yitSR data alonei- 

Equation ([3]) was fit to the data of Fig. [2] using the 
temperature dependence of the BCS order parameter 
for Ae(r) (Ref. 0), and varying A„ and the amplitude 
Ae(0) of Ae(T) for best fit. Figure [3] shows the depen- 
dence of Ae(0)/7p on solute concentration. The initial 
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FIG. 3. (Color online) Dependence of the rms width Ae (0) /7,i 
of the T = spontaneous field distribution on Ru concentra- 
tion X and La concentration y in Pr(0si_2,Rui)4Sbi2 and 
Pri_yLayOs4Sbi2. Solid lines: linear fits. Dashed curves: 
estimated experimental error. 



suppression is accurately linear for both solutes, with 
slopes -1.71(5) G (Ru doping) and -1.13(7) G (La dop- 
ing). For Ru doping the data are limited to a; < 0.3, and 
hence do not probe the crossover in penetration-depth 
behavior observed at higher concentrations.'^^ The dras- 
tic decrease in the specific heat jump at Tc observed in 
Pri_yLayOs4Sbi2 for y > 0.3 (Ref. HI) is not reflected in 
our data (Fig.[3|) [nor, for that matter, in Tc{y)]. 

The estimated error in Ae(0) diverges as Ae(0) — > 
(dashed curves in Fig. Thus data with experimen- 
tally attainable statistics (which are excellent, cf. Fig. [T]) 
cannot determine whether or not linearity is maintained 
at higher concentrations. Nevertheless the available data 
strongly suggest that broken TRS is suppressed for x > 
0.6 by Ru doping, but persists to high La concentrations. 

This behavior can be understood if the Cooper pair- 
ing mechanism in PrOs4Sbi2 is the exchange of itin- 
erant Pr'^"'" crystal-field excitations (excitons) 
Treatments of this interaction have concluded that 
it can lead to TRS-breaking superconductivity in 
PrOs4Sbisr^ and its alloys.— An alternative picture for 
Pr(Osi_a;Rua;)4Sbi2 explains the minimum in Tc as 
due to a crossing between the CEF splitting and the rat- 
tling energy of the Pr ion^^ without considering broken 



TRS. 

In the model of Koga, Matsumoto, and Shiba (KMS),— 
Ru and La doping of PrOs4Sbi2 affect the TRS-breaking 
excitonic pairing in different ways. The rapid increase 
with Ru doping of the splitting between the crystalline- 
electric-field (CEF) Pr''+ singlet ground state and mag- 
netic first excited state'^^ weakens the pairing interaction 
due to the excitonic mechanism, without pair breaking or 
other effects that would rapidly suppress Tc- PrRu4Sbi2 
is a conventional superconductor, and this weakening 
leads to a crossover or transition between TRS-breaking 
superconductivity in Os-rich alloys and conventional s- 
wave pairing at the Ru-rich end of the series. This is re- 
flected in the minimum in Tdx)^^ and, as reported here, 
in the vanishing of the broken TRS, both for x in the 
neighborhood of 0.6. 

In Pri_j,LayOs4Sbi2 La substitutes for Pr with little 
distortion of the lattice or the electronic structure J^i^^i^ 
La doping simply weakens the Pr-Pr intersite interac- 
tion, resulting in less excitonic dispersion; this reduces 
the pairing interaction^^ Our ZF-^SR data suggest that 
in this case TRS is broken across the alloy series, with 
an amplitude that vanishes only for large y. 

The temperature dependence of the exponential damp- 
ing rate A is given in Fig. 2] for the Ru-doped and La- 
doped alloys. As in PrOs4Sbi2,'' in the alloys there is 
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FIG. 4. (Color online) Temperature dependence of the ZF 
exponential damping rate A in (a) Pr(Osi-a;Rua;)4Sbi2, x = 
0, 0.1, 0.2, and 0.3, and (b) Pri_j;Laj;Os4Sbi2, y = 0, 0.2, 0.4, 
0.6. Arrows: Tc from bulk measurements. 



no evidence for an anomaly in A at Tc- Previous experi- 
ments"^^ showed that the dependence of the damping on 
longitudinal field is consistent with dynamic relaxation 
due to thermal fluctuations. Nuclear magnetism was sug- 
gested as the origin of these fluctuations for a number of 
reasons, among them the fact that electronic spin fluctu- 
ations would be strongly affected by superconductivity. 
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We consider for completeness a transition to an un- 
related weak-moment (< 10~^/xb) magnetic state at 
T = T^ng as an alternative explanation of B^. From the 
ZF-/iSR data the onset of Bg occurs at the superconduct- 
ing Tc, at least for low doping where it can be clearly seen 
(Fig. [2]). Although there are many cases of coexistence of 
magnetism and superconductivity in strongly-correlated 
electron systems, Tmag = Tc only at isolated points in the 
phase diagrams of these systems. A magnetic-transition 
scenario requires fine tuning to such a point. Even if this 
were the case in the end compound, doping would almost 
certainly change Tmag relative to and the fine tuning 
would be lost in the alloys. There is no evidence for this, 
and we conclude that a magnetic transition unrelated to 
superconductivity is unlikely. 

In Pr-based compounds the /i+ charge can affect the 
CEF splitting of Pr'^'^ near neighbors, which in turn 
modifies the local Pr'^^ susceptibility that is the major 
contribution to the muon Knight shift. '^^ This modifica- 
tion might also affect the superconductivity locally. The 
normal-state /z"*" Knight shift in PrOs4Sbi2 tracks the 
bulk susceptibility, however;^ suggesting that any such 
perturbation is small. 



We conclude that broken TRS in PrOs4Sbi2 is sup- 
pressed by both Ru and La doping, but differently for 
the two solutes. Ru doping appears to restore TRS for 
X > 0.6, near the minimum in Tc{x), whereas for La 
doping TRS breaking persists to y ^ 1, i.e., most or all 
Pr-doped LaOs4Sbi2 alloys exhibit broken TRS. These 
properties are consistent with the KMS picture^® for 
the CEF excitonic pairing mechanism and TRS-breaking 
superconductivity in PrOs4Sbi2-based alloys. Our re- 
sults motivate a quantitative treatment of broken TRS 
in Pr(Osi_a;Ruj;)4Sbi2 and Pri_j,LayOs4Sbi2. 
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